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ABSTRACT 

This paper presents a solely thermal flare, which we detected in the microwave range from 
the thermal gyro- and free-free emission it produced. An advantage of analyzing thermal gyro 
emission is its unique ability to precisely yield the magnetic field in the radiating volume. When 
combined with observationally-deduced plasma density and temperature, these magnetic field 
measurements offer a straightforward way of tracking evolution of the magnetic and thermal 
energies in the flare. For the event described here, the magnetic energy density in the radio- 
emitting volume declines over the flare rise phase, then stays roughly constant during the 
extended peak phase, but recovers to the original level over the decay phase. At the stage 
where the magnetic energy density decreases, the thermal energy density increases; however, 
this increase is insufficient, by roughly an order of magnitude, to compensate for the magnetic 
energy decrease. When the magnetic energy release is over, the source parameters come back 
to nearly their original values. We discuss possible scenarios to explain this behavior. 

Subject headings: Sun: flares—acceleration of particles—turbulence—diffusion—Sun: magnetic 
fields—Sun: radio radiation 


1. INTRODUCTION 

It is well established that to produce a flare requires 
magnetic free energy stored in the form of nonpo¬ 
tential coronal magnetic field (supported by coronal 
electric currents), which is somehow released and con¬ 
verted into other forms of energy—kinetic, thermal, 
and nonthermal. This process of magnetic energy re¬ 
lease is commonly referred to as magnetic reconnec¬ 
tion. 

Flaring energy release is necessarily an extended 
process. Indeed, to supply a large flare with an en¬ 
ergy of ^ 10^^ erg at the expense of magnetic energy 
dissipation, one needs a rather large volume. For ex¬ 
ample, ^ 10^^ erg is roughly the energy contained in a 
100 G magnetic field in a volume r\j 3-10^^ cm^. Since 
only free (nonpotential) magnetic energy is available 
for release, the energy release region must occupy an 
even larger volume. 

If the magnetic reconnection/energy release occurs 
in a given limited volume, it is reasonable to antici¬ 
pate that the magnetic field strength, energy density, 
and the total energy all decrease during the energy re¬ 
lease stage. Thus, to probe the very process of energy 
release and transformation to other forms, it would be 
highly desirable to track all the energy components in 
the flaring vol ume. This is. how ever, a highly com¬ 
plicated task (jEmslie et al.ll2QT^ . Indeed, although 
the thermal energy could often be reliably estimated 
with the soft X-ray (SXR) emission and some other 
means, the other energy components are very difficult 
to precisely measure at the flare energy-release site. 

The thermal flare energy can often be understood 
as the plasma response to accelerated particle input; 
the relationship commonly referred to as the Neu- 
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pert effect (jNenpertl 119681 ). However, there are cases 
(cold flares) when no detectable heating is observed 
in spite of highly efficient acceleration. In partic¬ 
ular, i n the 30 Jul 2002 col d flare we studied re¬ 
cently (jFleishman et al.ll2Qllf ) we found that almost 
all available electrons were accelerated. It was not 
possible, however, to construct the complete energy 
balance, because the fast electrons deposited their en¬ 
ergy into the chromosphere, as confirmed by footpoint 
hard X-ray (HXR) emission, although no significant 
evaporation occurred, and diagnostics of the chromo¬ 
spheric response could not be measured. From the gy- 
rosynchrotron (GS) spectrum we estimated the mag¬ 
netic field in the particle acceleration region, but ob¬ 
tained a magnetic field value that did not show any 
significant evolution. This can mean that, because 
the nonthermal electrons leave the acceleration re¬ 
gion relatively quickly, and thus do not contribute to 
the total pressure in the flaring volume, the magnetic 
field is supported by an external pressure at roughly 
its original level. 

On the other hand, in many events there is 
a prefiare phase th at is almost purely thermal 
(|Battaglia et al.ir2QQ9f ) , although sho wing a weak non¬ 
thermal component in many cases (jAsai et al.ll2QQ6l . 
l2QQ9l : lAltvntsev et al.l |2Q12[ ) . Since the nonthermal 
energy content is only minor over these preflare heat¬ 
ing episodes, the nonthermal particles alone cannot 
account for the heating. Thus, the plasma must be 
somehow heated directly by the magnetic energy re¬ 
lease, or indirectly, e.g., via dissipation of the turbu¬ 
lent motions and/or waves generated, perhaps, as a 
result of the magnetic energy release. But this implies 
that there can be cases when the nonthermal compo¬ 
nent is essentially nonexistent over the entire event 
duration, which would be manifested as a purely ther¬ 
mal event. 

Indeed, such cases have been observed in both 
microwave and X-ray domains. In particular. 
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iGarv fc Hurfordl (|1989[ ) reported a simple microwave 
burst showing the spectral shape consistent with be¬ 
ing produced by a purely thermal plasma (with a 
Maxwellian distribution without any nonthermal tail) 
over the entire evolution of the flare. From the se¬ 
quence of spectral fits they evaluated the character¬ 
istic magnetic field at the radio source to b e roughly 
consta nt at the level of ^ 770 G. Recently, iLiu et aP 
(l20T^ reported a thermal flare that occurred in a 
low-lying coronal loop. They derived evolution of the 
plasma temperature, but no direct data on the mag¬ 
netic field at the source was available due to lack of 
radio observations. 

Microwave data offer diagnostics of both magnetic 
field and thermal plasma. For a purely thermal flare 
the microwave emission is produced by a combination 
of the thermal GS plus free-free emissions. As an 
example. Figure [T] displays microwave spectra from 
a given volume uniformly filled with a hot plasma 
with various temperatures and having uniform mag¬ 
netic field. The exact spectra (thin lines) consist 
of a number of peaks corresponding to gyroharmon- 
ics (integer multiples of the gyrofrequency), whose 
flux level increases with frequency at low frequencies, 
where the emission at each gyroharmonic is optically 
thick and the corresponding brightness temperature 
is just the plasma kinetic temperature; thus the av¬ 
erage flux level rises as oc /^. Then, above a certain 
frequency, which increases as the temperature and/or 
magnetic field increases, the emission becomes opti¬ 
cally thin and so the flux level decreases quickly; thus, 
at even higher frequencies the flux level is determined 
by the optically thin free-free emission, which yields 
a flat spectrum at high frequencies. Isolated gyro- 
harmonic emis sion has been reported for non-flarin g 
active regions (iLang et al.lll98'^: iBogod et al.l|2QQQ[) . 
and from solar flares ( Benka fc HolmanI 119921) . but 
such reports are rare because the magnetic field is 
spatially nonuniform at the radio source. Even a 
relatively minor nonuniformity of the magnetic field 
will smooth the isolated harmonics, thus the thin line 
curves in Figure[T]are essentially unobservable in most 
cases. Instead, a continuum thermal spectrum is ob¬ 
served (thick lines in Figure [1]), with flux density oc 
at low frequencies, where the emission is optically 
thick, and a rapid (exponent-like) drop-off at higher 
frequencies, where the emission is optically thin; the 
falling part of the gyro-spectrum then gives way to a 
fiat free-free spectrum. 

The spectral peak of the thermal GS emission al¬ 
lows a precise measurement of the maximum mag¬ 
netic field at the source, while the level of the op¬ 
tically thick part of the spectrum depends only on 
the product of the plasma temperature and source 
area. Above the spectral peak the microwave spec¬ 
trum is fiat due to the optically thin free-free emission 
at the level defined by the plasma density and tem¬ 
perature. Thus, microwave observations of a purely 
thermal source offer a straightforward way of estimat¬ 
ing physical parameters needed to s tudy the energy 
balanc e and evolution in such events. lAltvntsev et al.l 
(|2Q12[ ) have demonstrated the ability of thermal GS 
emission to yield magnetic field diagnostics at the 
energy release site. However, they used microwave 



Fig. 1.— Model GS + free-free spectra from a uniform 
plasma volume computed with an exact expressions (specifi- 
cally, the hybrid fast GS code, IFleishman fc Kuznetsovll201CII . 
thin curves) and continuous approximation (thick curves). The 
adopted parameters are: area A = 1600"^, depth d = 10", 
plasma density no = 10^*^ cm“^, and viewing angle 9 = 75°. 
Values of the plasma temperature T and magnetic field B are 
shown in the panels. Note, that for a given magnetic field, 
the spectrum peak increases when the temperature increases. 
Thus, having a constant spectral peak while the plasma heating 
requires a corresponding decrease of the magnetic filed. 

data from a few discrete frequencies provided by the 
Nobeyama Radio Polarimeters (NoRP), RSTN, and 
Siberian Solar Radio Telescope (SSRT). As a result, 
the microwave spectral peak could not be precisely 
determined, which is needed for precise magnetic field 
diagnostics. 

Here we use spectrally resolved Owens Valley So¬ 
lar Array (OVSA) observations that provide us with 
much better measurements of the spectral shape suf¬ 
ficient to firmly constrain the source parameters such 
as temperature by area product and the magnetic 
field. We also use these exceptionally clean spectra as 
an opportunity to investigate the potential influence 
of a slight non-Maxwellian distribution, in the form 
of the kappa distribution. 

2. OBSERVATIONS 

The period of thermal emission that we study man¬ 
ifested during the decay phase of one of the largest 
flares of solar cycle 23, the GOES soft X-ray class 
X14.4 event of 2001 Apr 15 that peaked ^ 13:45 UT 
(Eigure [2^), and was spatially associated with the 
post-flare loops of that event. However, we believe 
the thermal event should more-properly be consid¬ 
ered as a ’precursor’ of a much smaller event that 
peaked ^ 17:00 UT (Eigure [2l3-c). Eigure [2]3 shows 
the Stokes I total power dynamic spectrum at 16:38- 
17:30 UT frorn the Owens Valley Solar Arra y (OVSA; 
iHurford et al. I [1983: IGarv fc Hurfordl 1 19941) . obtained 
from 40 frequencies roughly logarithmically spaced 
over the range 1-18 GHz. To aid in visibility of the 
emission, the corresponding 2.4 GHz Stokes I total 
power light curve is shown in Eigure [2t. No signifi¬ 
cant circular polarization was seen during the event. 
The spectral shape of the emission, as we will show 
shortly, is consistent with purely thermal gyrosyn- 
chrotron emission during the small, isolated peak to 
the left of the dashed line in Eigure [Sb-c (16:40-16:48 
UT), then it becomes increasingly nonthermal during 
the time to the right of the dashed line, leading up to 
the event peaking at 17:00 UT. Eigure [2]i-e show the 
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OVSA dynamic spectrum and 2.4 GHz light curve for 
the period of the thermal flare, from 16:38-16:48 UT. 
The narrow-band nature of the emission in Figure [2]i 
is striking. One sometimes sees highly nonthermal, 
coherent plas ma emission in th e decimeter range be¬ 
low 3 GHz (jNita et al.l l2QQ4l ). which OVSA’s fre¬ 
quency and time resolution may not resolve. How¬ 
ever, Figure [2]F shows much higher frequency- and 
time-resolution dat a from the Ondrejov observatory 
(|Jiricka et al.l 120011 ). on the same time scale as Fig¬ 
ure [2]i-e. The Ondrejov data conhrm the smooth con¬ 
tinuum nature of the emission, consistent with our 
interpretation of thermal gyrosynchrotron emission, 
and eliminates the possibility that the radio spectrum 
is affected by coherent plasma emission, which would 
be a signature of nonthermal particles. 

2.1. Context Observations 

The Tra nsition Region An d Goronal Explorer 
(TRAGE; iHandv et al.l 11999) and Solar Helio¬ 
spheric Observatory (SoH O; iDomingo et al.l 119951: 
iDelaboudiniere et al.l 1199'^ targeted this AR and re¬ 
vealed a system of postflare coronal loops from the 
XI4.4 event. The radio burst originated from the 
same area, but it is not possible to associate the 
radio source with a given single EUV loop, and in 
fact, as we will see, the radio source is far hotter 
than these EUV-emitting post-flare loops. Eigure[3t^- 
c show three of the TRAGE EUV wavelengths, 171 A 
(panel a), 195 A (panel b), and 284 A (panel c). The 
TRAGE images were greatly affected by high-energy 
particle noise from the earlier XI4.4 flare, hence we 
have processed the images with the eit.despeckle rou¬ 
tine from the Solar Software (SSW) IDL package. 
The images show a system of post-flare loops with 
bright tops centrally located over the active region, 
with more diffuse loops to both the north and south. 
Both footpoints of the post-flare loops appear to be 
on the disk, just inside the limb. 

Hydrogen-alpha images from Big Bear Solar Ob¬ 
servatory (BBSO) are available at 16:00 and 17:00 
UT. We show in Eigure [3li the relevant portion of 
the 16:00 UT Ha image, showing an earlier stage of 
post-flare loops. The BBSO Ha image at 17:00 UT 
shows that the loops with the bright-looptops have 
migrated southward by then. 

2.2. OVSA Imaging and Relative Visibility 

As shown in Eigure [2l the OVSA radio burst has 
signihcant flux density only in the range 1.2-3.2 GHz. 
Eigure [3] displays a summary of the OVSA imag¬ 
ing at 1.2-3.2 GHz at the peak time of the flare, 
^16:42:35 UT, overlaid on the four panels mentioned 
above. Panels a-c show contours at 20, 40, 60 and 
80% of the peak intensity for 1.2 GHz (panel a), 
1.4 GHz (panel b), and 1.8 GHz (panel c). The 
dashed ellipses in each panel show the sythesized half¬ 
power beam size. Eigure [3li shows the relationships 
of source position and size for the range of frequen¬ 
cies of the burst, as 50% contours at 1.2 GHz (red), 
1.4 GHz (orange), 1.8 GHz (yellow), 2.4 GHz (green) 
and 3.2 GHz (blue); see the on-line movie for greater 
detail. The half-power synthesized beam (dashed el¬ 
lipse) shown in this panel is for the frequency 2.4 GHz. 
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Fig. 2.— The thermal flare event in the context of other 
radio activity, (a) The GOES (1-8 A) soft X-ray light curve 
spanning 12 hours (1200-2400 UT) on 2001 Apr 15, showing 
the large (X14.4) flare peaking around 13:50 UT and its long 
decay. The vertical dashed lines show the timerange shown in 
b-c. (b) The OVSA dynamic spectrum from 1-18 GHz from 
1640-1730 UT, showing the narrowband thermal flare (left of 
the white dashed line) and an increasingly nonthermal preflare 
phase leading to a larger event peaking at 1700 UT. (c) The 
2.4 GHz lightcurve corresponding to b. To the left of the ver¬ 
tical dashed line is the timerange of the thermal flare shown 
in d-e. (d) The OVSA dynamic spectrum from 1-18 GHz from 
1638-1648 UT, showing the narrow-band thermal flare, visible 
only in the range 1.0-3.2 GHz. (e) The 2.4 GHz lightcurve 
corresponding to d. (f) The 1-2 GHz dynamic spectrum from 
Ondrejov on the same timescale as in d-e. The horizontal and 
vertical lines are instrumental artifacts. 


We have studied the time evolution of the burst at 
each frequency, and find no significant change of po¬ 
sition; see the on-line movie for greater detail. The 
sources at each frequency rise and fall in intensity 
from a single location over the 4-5 minute duration of 
the burst. 

The radio sources are spatially associated with the 
centrally-located post-flare loops-those with bright 
EUV emission at the loop-tops. As we will see 
shortly, the radio spectrum of the burst is con¬ 
sistent with purely thermal gyrosynchrotron emis¬ 
sion from a homogeneous source. However, such a 
source should have a fixed location and frequency- 
independent source size, seemingly at odds with the 
frequency-dependent source sizes and locations shown 
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TABLE 1 

Table 1: Comparison of deconvolved source size for 

OBSERVATIONS AND SIMULATION 
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Fig. 3.— Radio images and other context imaging data, (a) 
1.2 GHz OVSA (contours) on top of TRACE 171 A (image), 
(b) 1.4 GHz OVSA (contours) on top of TRACE 195 A (im¬ 
age), (c) 1.8 GHz OVSA (contours) on top of TRACE 284 A 
(image), and (d) OVSA 50% contour at each of 1.2 (red), 1.4 
(orange), 1.8 (yellow), 2.4 (green) and 3.2 GHz (blue), on top 
of BBSO Hq; (image), showing a slight dispersion in position 
with frequency. See on-line movies showing the source varia¬ 
tion either with frequency or with time. 

in Figure [3li. Of course, much of the source variation 
could be due to convolution of the source with the 
frequency-dependent synthesized beam (point spread 
function). To check this, we use d the s imple analyt¬ 
ical deconvolution approach of IWildl (jl97Qf ) , which 
allowed us to use the known 2d-gaussian parame¬ 
ters of the synthesized beam and those of the ob¬ 
served source to derive the parameters of the decon¬ 
volved source (semi-major axes a and 6, and orienta¬ 
tion), from which the deconvolved source area irah 
could be determined. The results of this analysis 
still yield a frequency-dependent deconvolved source 
size (’’observed” column of Table [T]). However, the 
OVSA interferometric array comprises only 6 anten¬ 
nas, one of which suffered from instrumental effects 
and was flagged from the data. Instantaneous (snap¬ 
shot) imaging with such a small number of anten¬ 
nas (5) provides only a small number of constraints 
on source parameters, which may be insufficient for 
an unambiguous determination of deconvolved source 
size. To check this possibility, we used the uvmo del 
routine from the Miriad package (jSanlt et al.ll2Qlll ) to 
create a model visibility database using a 2-d gaussian 
source with fixed source size vs. frequency (we chose 
the deconvolved 1.8 GHz source area nah = 4113 sq. 
arcsec, with a = 39.1" and b = 33.5"). We then 
sampled the model visibilities with the OVSA an¬ 
tenna spacings. This model database was then used 
to create simulated OVSA images, which were decon¬ 
volved in the same way as the original OVSA images. 
Table [T] shows the comparison of the original decon¬ 
volved source area vs. frequency with that of the 
simulated images. It is clear that the same decrease 
of source area with frequency is seen in both cases, 
despite the fact that the simulation source model has 
a constant area. This implies that determining the 


Preq. 

(GHz) 

Size Observed 
(arcsec^) 

Size Simulated 
(arcsec^) 

(j) Observed 
(degrees) 

1.2 

12742 

11097 

346 

1.4 

7241 

7633 

343 

1.6 

5252 

5013 

345 

1.8 

4113 

3643 

341 

2.0 

3552 

2995 

329 

2.4 

2301 

2158 

348 

2.6 

1736 

1780 

355 

2.8 

1205 

1568 

350 


source area from the 6-antenna OVRO imaging ob¬ 
servations is not reliable, but also that the observa¬ 
tions do not rule out a roughly constant source size 
with frequency, as expected for a therm al so urce. We 
will come back to this point in sections 12.31 & IH 

Although we have found that the deconvolved 
source size is not reliable, it is interesting that the 
orientation of the convolved (observed) sources in 
Figure [3] is slig htly d i fferent from the orientation of 
the beam. The IWildl (jl97Qf ) deconvolution yields an 
orientation 0 for the deconvolved source of around 
345 deg (fourth column of Tabled]), which is nearly 
perpendicular to the limb, i.e., in the same direc¬ 
tion as the slight dispersion in centroid height with 
frequency shown in Figure [3]i. This may suggest a 
vertically-oriented hot loop with a slight frequency 
asymmetry, i.e. low frequencies tending to come from 
near the top and high frequencies from near the bot¬ 
tom. However, as we shall see later, any such inho¬ 
mogeneity must be slight enough not to broaden the 
total power radio spectrum significantly. Stated an¬ 
other way, although the imaging data imply that the 
sources at different frequencies are displaced relative 
to each other, they must all have roughly the same 
effective areas. 

Given the deconvolved source sizes are not reliable, 
we can study the frequency dependence of the relative 
visibilities to determine the source size assuming that 
there is a single source at all frequencies, as in an other 
similar case reported by iGarv Sz HnrfordI (jl989f ) . Re¬ 
call that the relative visibility is the ratio of baseline 
amplitude to the square-root of the product of the to¬ 
tal power on the two antennas corresponding to that 
baseline. It has the advantage that it is independent 
of calibration. Wha t one expects to see fo r a source 
of constant size (see IGarv fc HnrfordI 1 1989f ) is a con¬ 
stant negative slope vs. frequency^, with the slope 
proportional to the size of the source. Gontrary to 
expectations, however, this procedure yields on all 
baselines a steadily rising relative visibility to 3 GHz 
(points above this are just noise, since the spectrum 
dies above this). If taken at face value, this implies a 
spatially-complex source. The relative visibility spec¬ 
tral shape is remarkably constant in time throughout 
the burst, so there is no indication that the source is 
changing size or shape with time. 

2.3. Understanding the spatial relationships 

Putting all of the evidence together, we get a pic¬ 
ture of a spectral shape and temporal behavior con- 
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IRCP+LCP][s1ul ® 1.2 GHz 


[RCP+LCPl[sfu] @ 1.6 GHz 


[RCP+LCP][s1u] ® 2.2 GHz 


[RCP+LCP]Is1u] ® 3.6 GHz 



Fig. 4.— A model thermal loop constructed to prove the feasibility of the ’gyro-stripe’ model concept. We used a nonlinear 
force-free extrapolated magnetic data cube obtained for AR 121 58 based on the H MI vector magnetogram taken close to the disk 
center on September 10, 2014. Using the GX Simulator tool (INita et al.l 12014( 1 we selected a magnetic flux tube, rotated the 
data cube to the western limb, and restricted the held of view to observe the emission from a fraction of the loop seen in the 
upper left panel. Then, the e mission was computed using the ’hybrid’ (containing the gyroharmonics) version of the fast code 
(IFleishman &: Kuznetso^1120101) to obtain the spectrum (upper right) and the radio images at different frequencies. The set of 
images (see on-line movie for better impression) proves that the sources at different frequencies can be (i) complex, i.e., consist of 
a few gyro-stripes embedded in a weaker and more diffuse continuum, (ii) displaced relative to each other, while (iii) have roughly 
same area and (iv) produce an integrated spectrum (upper right) that nicely fits the observed spectrum, thus, emulating a uniform 
thermal source. Note that in a general case the model integrated spectrum will be broader than the observed narrow quasi-uniform 
thermal spectrum. 


sistent with a relatively compact source of uniform 
area, but with contrary spatial information that in¬ 
stead indicates a source with significant spatial com¬ 
plexity. We have insufficient information to be defini¬ 
tive as to how to resolve this quandary, but one pos¬ 
sibility is that we are observing the manifestation of 
discrete gyro-layers expected for non-uniform sources, 
coming from the harmonic spectral features shown in 
Figure [H Such feature s have been suggested before 
(jKnzne^ov et al.ll2Qll[ ). as shown in the movie ac¬ 
companying the online version of this paper, some 
frames of which are shown in Figure IH Under cer¬ 
tain special conditions (which could account for the 
rarity of such events), the total emitting area of such 
a non-uniform source, made up of up to 5-7 gyro- 
layers, could be remarkably constant with frequency 
over some restricted range, such as the 1.2-3.2 GHz 
range of the present event, but the centroid of the 
emission at different frequencies can be displaced fol¬ 
lowing the apparent spatial displacement of the gyro 
stripes as a function of frequency. Generally, con¬ 
tributions from two different legs of the loop can 
broaden the spectrum. However, for a loop observed 
’edge-on’ at the limb—the geometry supported by the 


recovered source position angle and the source dis¬ 
placement with frequency—the farthest leg can be 
occulted by the closest one; so only the leg closest to 
the observer contributes to the spectrum. We return 
to this point later. Attempting detailed modeling of 
this, however, is not well constrained by the data and 
so beyond the scope of this paper, but the proof-of- 
concept test in Figure |4] suggests that it is feasible. 

3. RADIO SPEGTRAL FIT 
3.1. Fitting framework 

In most of the time frames (see symbols with er¬ 
ror bars in Figure [5j) the microwave spectrum of the 
burst rises with frequency oc at low frequencies 
until terminated by a very sharp, almost exponential 
cutoff, which then gives way to a roughly fiat flux 
level at the correspondingly higher frequencies. Such 
a spectrum can be fit perfectly by a combination of 
gyrosynchrotron and free-free emissions from a purely 
thermal source with quasi-uniform area as a function 
of frequency, where the gyro emission is responsible 
for the low-frequency, relatively narrowband spectral 
peak, while the optically thin free-free emission is re¬ 
sponsible for the higher-frequency broadband back- 
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ground. 

Within this spectral model, the low-frequency re¬ 
gion, oc /^, represents the optically thick thermal 
gyro emission, thus, the flux level is defined solely 
by the product of the source area A and the plasma 
temperature T, F oc A • T. In particular, know¬ 
ing the plasma temperature, e.g., f rom X-ray data 
allow s estimating the source area (|Altvntsev et al.l 
|2Q12[ ). while having information on the source size 
allows reading the plasma temperatur e from the spec¬ 
trum directly (|Garv fc Hur for dill9891) . There are no 
X-ray data available to derive the plasma temper¬ 
ature for this event, so we must rely on the other 
constrai nts for the source size. As has been said in 
Section 12.2[ the estimate of the source area is not 
reliable. It is important, however, that the imag¬ 
ing data do not suggest any noticeable change in 
the source size or position with time, while at the 
same time the narrow spectral shape suggests that 
the source area cannot vary too much with frequency, 
as any appreciable variation will broaden the radio 
spectrum. In addition, analysis of the relative vis¬ 
ibilities suggests a noticeable spatial complexity of 
the radio source, which, all together, drove us to a 
’gyro-stripe’ model envisioned in Sec 12.31 thus, im¬ 
plying the effective source area to comprise only a 
fraction of the deconvolved source area. Accordingly, 
we loosely adopt an area of Am = 1600"^, which is 
roughly 40% of the deconvolved source area found at 
1.8 GHz, Aq = 7r(39.1" x 33.5") = 4113"^. One has to 
keep in mind that any temperatures derived from this 
assumption should be scaled by A^/A, where A is 
the true source area. The source depth is adopted to 
be 20", which is an appropriately scaled deconvolved 
source width (the exact value of the source depth is 
unimportant; for example, very similar results would 
be obtained for 10" depth, see Figured]). 

A striking feature of this radio burst is nearly per¬ 
fect constancy of the spectral peak. Indeed, for a 
thermal source with constant size, the increase of 
the flux level can only be provided by an increase 
in plasma temperature. In a source with a constant 
magnetic field, this plasma heating must necessarily 
result in the corresponding increase of the spectrum 
peak, see Figure dl Thus, the aforementioned con¬ 
stancy of the spectral peak at the radio emitting vol¬ 
ume implies a corresponding decrease of the magnetic 
field over the course of the plasma heating. 

3.2. Fit to Maxwellian Plasma 

To perform th e sequential s pectr al _ fit 

(iBastian et al.l I2QQ7I: iFleishman et al.l [20091 120131: 
IGarv et a.lJ l2nTl we apply a homogeneous source 
model with constant depth 20" and area 1600"^ 
and three unknown free parameters H(t), T(t), 
and the magnetic field, plasma temperature, 

and density respectively. Selection of the objective 
function to be forward fit requires some further 
discussion in this case. Indeed, the gyro emission 
from a truly homogeneous thermal source (with a 
unique single magnetic field value), called gyroreso- 
nance (GR) emission, would represent a number of 
distinct narrow peaks occurring around a few small 



10 10 10 10 
Frequency, GHz Frequency, GHz Frequency, GHz Frequency, GHz 


Fig. 5. — Examples of the spectrum fits. See the on-line set 
for the entire flare duration. Only two significant digits of k 
are kept. 

integer multiples of the gyrofrequency (gyroharmon- 
ics). Figure dl although a realistic variation of the 
magnetic field along the line of sight smooths these 
peaks out to form a continuum spectrum. Formation 
of such a continuous spectrum is accounted for by 
the theory of GR emission, which properly takes 
into account the ma gnetic field variati on along 
the line of sight (e.g.. I^eleznvak^ ll97Qf ) and has 
been widely applied to the s ola r active regions 


(lAlissandrakis & Kundu 

119841: 

Akhmedov et al. 

1986 


IGarv & Hu 

jford 

1994 2004 

: iPeterova et al. 

2006 

iLed 120071: 

BoKod 

& Yasnov 

20091: iTun et al. 

2011 

Ixita et al.l 

2®^ 

Kaltman et al.l 1201211. For 


the purpose of the practical forward fit there is no 
need to include the magnetic field nonuniformity 
explicitly. Instead, we can use the continuous version 
of the fast GS code ([Fleishman fc Kuznetsov! 1^01 Q[ ) 
from a homogeneous source, which automatically 
performs the required averaging. We checked via 
modeling that the results of these two approaches 
agree very well to each other. 

We performed the fit of the observed sequence of 
the microwave spectra to the continuous fast GS 
code with the Maxwellian distribution of the ra¬ 
diating electrons using the simplex method with 
shaking, which proved to work remarkably well for 
the microwave sp ectra ([Fleishman et al.l I2QQ9L l2Q13l : 
IGarv et~aIll2Q13[ ). The weights of the individual data 
points were calculated taking into account the mis¬ 
matches between different OVSA antenna measure¬ 
ments at the given frequency and time frame. As 
has been said, there are three free parameters of the 
model, H(t), T(t), and whose evolution we re¬ 

covered from the sequential fit. The fit is successful 
throughout the entire burst duration; examples of the 
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fit are given in Fig. [5] by black solid curves. 

3.3. Fit to a Plasma with Kappa Distribution 

Interestingly, with this set of the microwave spec¬ 
tra we can address the question of how precisely 
the Maxwellian distribution is maintained in the ra¬ 
dio source. To do so we perform a similar sequen¬ 
tial fit b ut in this case assuming the kappa dis- 


tribntion (iVasvlinnasI!19681: lOwocki & Scndderl 

1983; 

Maksimovic et ah! 

ll997l:lLivadiotis & McComasI 

2009; 

Pierrard & Lazarll 

1201011 of the olasma instead of the 


Maxwellian one. The kappa distribution (of the sec¬ 
ond kind in the notation of iLivadiotis &: McComasI 
[200^ contains one more free parameter, the in¬ 
dex compared with the Maxwellian distribution, 
(for the microwa ve emission, produced by these 
distributions, see iChinderi fc Chinderi Dragg I2QQ4 
iFleishman fc KnznetsovI I2Q1QL |2Q14 and references 
therein). In a spatially uniform source, the value 
of the index k can be used as a measure of the de¬ 
viation of the actual plasma distribution from the 
Maxwellian one in such a way that the kappa distri¬ 
bution with large n is very similar to the Maxwellian 
one, while that with s mall n is marked ly different 
from it. For example, iBian et ^ (j2Ql4) have pro¬ 
posed that the kappa distribution can represent a so¬ 
lution of a transport equation for nonthe rmal elec¬ 
trons in a dense acceleration region, while lOka et all 
found that the X-ray spectra from above-the- 
loop-top (Masuda) sources can be well fit with the 
kappa distribution with indices in the range 4-14, 
alt hough even smaller n i ndices have been reported 
bv iKasparova fc KarlickA^ (j2QQ9f ) for the coronal loop- 
top sources. On the other hand, in a microwave total 
power spectrum, spatial nonuniformity can mimic the 
effect of finite kappa, i.e. the spectrum of a thermal 
source will be broadened by the nonuniformity com¬ 
pared with the thermal spectral shape discussed in 
section [321 We checked, via similar modeling to that 
presented in Fig.jH that the kappa distribution can of¬ 
ten provide an appropriate fit to the total power spec¬ 
trum from a nonuniform source, even though filled 
with a purely Maxwellian thermal plasma. Having in 
mind both of these possibilities, we apply the kappa- 
distribution fit phenomenologically, without immedi¬ 
ate implication for the non-Maxwellian nature of the 
plasma or the source nonuniformity. 

The fit examples with the kappa distribution are 
given by the solid green lines in Fig. [5] along with the 
n indices obtained for each time frame. Note, that 
for /i: > 15 — 20 the microwave spectra produced by 
either Maxwellian or kappa distribution are visually 
indistinguishable in Fig. [5] from each other—the cor¬ 
responding kappa area is marked yellow in the kappa- 
index panel of Fig. |6l stated another way the electron 
energy distribution in the radio source is compara¬ 
bly consistent with both the Maxwellian distribution 
and the kappa distribution with a reasonably large in¬ 
dex over most of the flare duration. However, at the 
late decay phase the fit with kappa distribution (with 
progressively smaller n indices) is formally preferable 
as it yields a significantly smaller measure. Nev¬ 
ertheless, we believe that this does not necessarily 
indicate a corresponding change of the (local) ther¬ 


mal electron energy distribution, but rather implies 
a progressively enhanced role of the source nonuni¬ 
formity due presumably to differing cooling speeds 
of the adjacent fluxtubes forming the radio source. 
This further implies that the kappa distribution in¬ 
stead of the Maxwellian one can successfully be used 
for analysis of the thermal flares with significant spa¬ 
tial nonuniformity. 

3.4. Fit Results 

Figure [6] displays evolution of the fit parameters 
H(t), T(t), and nth{t) (and the n index in case 
of kappa distribution) and also the energy densities 
along with a few radio light curves shown for the ref¬ 
erence. The derived plasma temperature, which is 
almost identical for both Maxwellian and kappa fits, 
displays clear rise, extended peak, and decay phases 
indicative of the plasma heating, saturation, and cool¬ 
ing. The plasma density displays a similar rise-peak- 
decay pattern although with noticeably larger scatter 
of the individual data points; the origin of this scat¬ 
ter is apparent: the number density estimate comes 
from the flat level of the microwave spectrum at high 
frequencies, where it is measured with large uncer¬ 
tainties. In contrast to these plasma parameters, the 
magnetic held displays remarkably different behav¬ 
ior: it goes down, stays relatively small, and then 
increases back to almost original values in case of 
Maxwellian fit, but in case of the kappa fit it remains 
at a relatively low level up to the late decay phase. 

Interestingly, to obtain a good fit a significant vari¬ 
ation of the magnetic held at the source is unavoid¬ 
able. It is tempting to conclude that we detect here 
the very process of the magnetic to thermal energy 
conversion; perhaps, due to the magnetic reconnec¬ 
tion process. However, such a straightforward con¬ 
clusion looks unwarranted: although the thermal en¬ 
ergy density does go up as the magnetic one declines 
at the flare rise phase, this is insufficient to make up 
for the magnetic energy density deficit. Furthermore, 
the very process of the magnetic to thermal energy 
conversion would imply the magnetic energy decrease 
when the plasma temperature keeps roughly constant 
during the extended peak of the burst to make up for 
the significant conductive losses (see next Section); 
however, the sum of the magnetic and thermal en¬ 
ergy keeps remarkably constant there, which requires 
another energy source to compensate for this loss. 

A more likely scenario comes from the loop ’self- 
occultation’ model envisioned in Sec ESI Let us con¬ 
sider a nested loop system, one above the other (with 
the outer ones having a progressively smaller mag¬ 
netic held), seen edge-on. Initially, when the tem¬ 
perature is small (coronal), the outer loops are trans¬ 
parent for the radio emission at 1.2-2.6 GHz (the ob¬ 
served optically thick part of the spectrum), so we see 
emission from the most inner loop with the largest 
magnetic held (^ 200 — 250 G). Then the tempera¬ 
ture in the outer loops start to grow and they become 
opaque sequentially—the larger the magnetic held in 
the given loop the sooner. This way the inner loops 
become occulted by the outer loop sequentially, so 
we observe emission from the most outer opaque loop; 
the hotter the loop the smaller the magnetic held that 
































Fig. 6. — Radio source parameters as derived from the OVSA 
spectral fit for three parameters of the thermal GS + free-free 
source. Two top panels show the radio light curves at a few 
frequencies indicated in the plots. The third right panel dis¬ 
plays the index of kappa-distribution derived from the fit. The 
yellow area indicates the /c values that give a distribution in¬ 
distinguishable from the Maxwellian one; the k, uncertainties 
are of the order of the k, values or larger in this area. Below the 
yellow area the formal uncertainties of the k index are about 
10% and are within the symbol size. The lower left panel shows 
evolution of the energy densities: magnetic (circles) and ther¬ 
mal (plus signs) energy densities as derived from Maxwellian 
fit and their sum (asterisks) as well as the sum derived from the 
kappa-distribution fit (green asterisks). The last panel shows 
the evolution for these two fits. 


makes it opaque at the given frequency range; this is 
why the fit gives us the proportionally smaller mag¬ 
netic field for the progressively larger plasma temper¬ 
ature. Apparently, for this model to work we have to 
adopt that the outer loops are heated progressively 
stronger than the inner loops. This model is capa¬ 
ble of explaining both constancy of the spectral peak 
during the heating-cooling evolution and constancy of 
the source location and area in time as observed. Al¬ 
though this ’fine-tuning’ of the proposed model may 
look artificial, we can argue that having a truly ex¬ 
ceptional event does require some really seldom com¬ 
bination of physical parameters and/or geometry. 

3.5. A Consistency Check 

The plasma cooling at the burst decay phase must 
be driven by either radiative or conductive cooling, 
whose characteristic time scales are strictly specified 
by the source sizes, temperature, and de nsity defined 
abov e. The radiative cooling time (e.g.. lAschwandenI 
l2QQ5f ) for such a hot plasma (T ^ 40 MK) is much 
longer than the observed one, r ^ 30 s; thus the 
conductive cooling is likely to dominate. The con¬ 
ductive co oling time, Trond ^ L^pCpIiKsT^/^) (Eq 
(4.3.10) in lAschwandenI [2005f ). where L is the loop 
length, p = UerUp is the plasma mass density, for 
the fully ionized hydrogen ideal gas Cpp = 
where Cp = 2^kB/[{'y — l)^p] erg/(g-K) is the spe¬ 


cific heat, 7 = 5/3, and ks ~ 9.2 • 10 ^ erg cm ^ 
s“^ is the Spitzer conductivity coefficient (see , 

e.g.. lAschwandenir2005l: iFleishman fc Toptvginll201^ . 
Solving for T and substituting known constants, we 
find 



which, for the observed source length L r\j 6 • 10^ cm 
and decay time scale r ^ 30 s and the plasma den¬ 
sity Ue ^ 10^^ cm“^ derived from the fit, yields the 
plasma temperature T ^ 40 MK in full agreement 
with the numbers derived from the fit. 

4. DISCUSSION 

The ability to precisely derive two components of 
the flare energy density—the magnetic and thermal 
ones—in the event under study (where the third gen¬ 
erally important energy component, the nonthermal 
one, is undetectable and so negligible) raises a ques¬ 
tion of how the energy is being transformed from one 
form to another and what are the means of the energy 
dissipation, escape, and build up. 

There are two possible scenarios of the event: (1) 
there is a single (isolated) volume where all these en¬ 
ergy transformations occur, in particular, where the 
magnetic energy is dissipated and then recovers back; 
or (2) the apparent source of the radio brightness 
moves in space; in particular, different magnetic field 
values obtained at different times pertain to physi¬ 
cally different spatial locations (layers) rather than 
reflecting the time evolution of the magnetic field at 
a given layer. 

4.1. Scenario One: A Single Isolated Source 

In case of an isolated source, there must be a bal¬ 
ance between its energy components, two of which are 
detected. The energy balance can be written in the 
form 

( 2 ) 

where u{t) is the fluid velocity, the integration is per¬ 
formed over the source volume that can change over 
the course of the flare in a general case, but stays 
constant (at least, the source area) in our event as 
suggested by OVSA imaging, and El is the total en¬ 
ergy loss rate, which includes the radiative and con¬ 
ductive losses; the factor of 3 in the third term is 
written in assumption of the thermal energy equipar- 
t it ion between electrons and ions. It is easy to esti¬ 
mate that both conductive and radiative losses can 
only remove a fraction of the thermal energy, which 
is itself small compared with the total energy domi¬ 
nated by the magnetic field. This means that with a 
good accuracy we can adopt 
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/ (^ + dv = 

( 3 ) 

Now, for a source with a constant volume we must 
require that the sum of the mean energy densities is 
conserved: 


BHt) 

Stt 


- ^—- + 3nth{t)kBT{t) 

« 2.3 • 10^ erg cm“^, 


r2 

Stt 

( 4 ) 


where Bq is the initial magnetic field value at the 
source. It is easy to estimate that the thermal energy 
density 1.6 • 10^ erg cm“^ at the peak phase) is 
a minor fraction (^7%) of the total energy density 
and so can be discarded in the zeroth approximation. 
Then, we note that the magnetic field at the peak 
phase of the burst is roughly two times smaller than 
the initial magnetic field value, see Figure [6l thus, the 
magnetic energy density (6.7 • 10^ erg cm“^) contains 
only about 30% of the original value so the remaining 
~70% must be ascribed to the kinetic energy density, 
pv? {2 ^ 1.6 • 10^ erg cm“^, within this model. Such 
an energy partition cannot be supplied, for example, 
by an ensemble of MHD waves: in the Alfven wave 
the kinetic and magnetic energy densities equal to 
each other; thus, together with the mean magnetic 
field component, the magnetic energy density must 
dominate in case of the MHD waves. The dominance 
of the kinetic energy would imply that the plasma 
motions are supersonic and super-Alfvenic, i.e., the 
shock waves. Although not firmly excluded, having 
such a strong kinetic velocity field in the entire vol¬ 
ume seems unlikely for many reasons. One of them is 
a clear lack of the radio image evolution, which would 
be expected in the presence of multiple shock waves in 
the source volume. In addition, such a strong kinetic 
velocity f ield would efficiently accel erate the charged 
particles (|Bvkov fc FleishmaniDQQQD , which is not ob¬ 
served either. 

Finally, there is an option of changing volume (al¬ 
though not supported by observations), when the ob¬ 
served decrease of the magnetic energy density is com¬ 
pensated by the corresponding increase of the flare 
volume, i.e., the volume expansion. The volume ex¬ 
pansion requires a corresponding net driving force. 
Conversion of the magnetic to thermal energy cannot 
produce this force since this conversion results in a 
decrease of the total pressure provided that the total 
energy is conserved. Thus, the volume expansion nec¬ 
essarily requires a decrease of the magnetic tensions, 
which can only be achieved through the magnetic re¬ 
connection process. This cannot be understood with¬ 
out considering a bigger volume than the immediate 
radio source. 


4.2. Scenario Two: an Open Source 

If we allow the source of the apparent radio bright¬ 
ness to move along the line of sight (recall, no dis¬ 
placement in the transverse direction is supported 
by observations) to illuminate different spatial lay¬ 
ers with respectively various magnetic fields detected 


at the different time frames, we have to consider the 
energy balance in a bigger volume in which, however, 
we do not have direct information about spatial distri¬ 
bution of the magnetic field and thermal plasma. In¬ 
stead, we consider three different stages of the burst— 
rise, peak, and decay. 

During the rise phase we clearly detect the pro¬ 
cess of the plasma heating during which the hotter 
the plasma the smaller the magnetic field at the ra¬ 
dio source. We have already noted difficulties with 
a model with magnetic dissipation at a given spa¬ 
tial location to account for this detected decrease of 
the magnetic field. Thus, we have to conclude that 
the process of plasma heating moves up from some 
low layer to higher coronal layers having correspond¬ 
ingly smaller magnetic field as observed. In our ’self- 
occulting’ source model the outer layers are originally 
cool and, so, optically thin. As soon as the outer 
layers with smaller magnetic field are heated they be¬ 
come opaque for the gyro process; thus, the lower lay¬ 
ers with the stronger magnetic field become invisible. 
This explains the observed evolution of the magnetic 
field. Apparently, the heating rate is stronger than 
the cooling rate at this rise phase. 

The extended peak phase of roughly flat light ra¬ 
dio curves deserves a close attention. Here, the ther¬ 
mal energy density, although small in the total en¬ 
ergy budget, is not small any longer compared with 
the magnetic energy density at the peak flare phase. 
Indeed, a constant temperature level at about 40 MK 
implies a relatively high heating rate to compensate 
rapid conductive losses with Tcond ^ 30 s. Given that 
the thermal energy density remains roughly constant 
at the level of ^ 160 erg cm“^ over about two min¬ 
utes, the source volume must have been supplied with 
^ 500 erg cm“^ over those two minutes. This value 
is comparable to the total magnetic energy density 
at the radio source during this extended peak phase. 
Nevertheless, the magnetic field keeps roughly con¬ 
stant during this phase, which unavoidably requires 
some energy input from an external energy source. 
Even more remarkable, the sum of the magnetic and 
thermal energies keeps precisely constant during the 
peak phase. 

When this energy source is switched off or ex¬ 
hausted, the conductive plasma cooling starts to con¬ 
trol the further evolution. The higher the tempera¬ 
ture the faster the cooling; thus, the effective radi¬ 
ating layer moves back down and the magnetic field 
increases again almost to the initial level. 

4.3. Energy Distribution of the Flare Plasma 

One more fundamental question that can be con¬ 
clusively addressed with this particular event, whose 
spectral evolution is consistent with a purely thermal 
process of the heating and cooling of the Maxwellian 
plasma, is what is the accuracy with which the 
Maxwellian distribution is maintained during the 
flare, or, equivalently, what deviations from the 
Maxwellian distribution are consistent with the data. 
To quantify these deviations we have considered, in 
Section [3l two spectral fits of the microwave spec¬ 
trum sequence—one with the Maxwellian distribution 
of the plasma electrons, and the other one with the 
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kappa distribution. A convenience of the kappa distri¬ 
bution for the purpose of characterization of presum¬ 
ably small deviations of the actual distribution from 
the Maxwellian one is that it depends on the n index 
in such a way that the kappa distribution is equiva¬ 
lent to the Maxwellian one for ^ oc, but develops 
a more and more pronounced power-law tail for pro¬ 
gressively smaller k. index. Other parameters of the 
kappa distribution are the same as for the Maxwellian 
one; it is important that the temperature T of the 
kappa distribution has the same physical meaning 
(the second moment of the distribution characterizing 
the mean energy of the electrons) as the temperature 
T of the Maxwellian distribution. 

We note that the measure^ for both fits are 
comparable to each other over the rise and peak 
phases of the burst. This implies that the data are 
comparably well fit by both these distributions. In 
the case of kappa-distribution fit, the k index is rel¬ 
atively large here, k > 100, with typical uncertain¬ 
ties in K larger than 100%. Such large values of k 
result in visually indistinguishable microwave spec¬ 
tra produced by either Maxwellian or kappa distribu¬ 
tion. One can conclude that the data are consistent 
with the kappa distribution with a large index; but 
a comparably meaningful conclusion is that there is 
no objective evidence in favor of any deviation from 
the Maxwellian distribution over the rise and peak 
phases. Importantly, the recovered physical param¬ 
eters of the source, 5(t), T(t), and are all 

consistent for both fits over these phases. 

The decay phase is apparently different from the 
rise and peak phases. Indeed, the is clearly 
smaller for the kappa distribution than for the 
Maxwellian one, which formally favors the kappa over 
the Maxwellian distribution. The hc index decreases 
progressively at the decay phase from 8 to 2, with 
the uncertainty about 10%, the magnetic field de¬ 
rived from the kappa fit is noticeably smaller than 
that derived from the Maxwellian fit; however, the 
temperatures still agree very well with each other. 

Thus, the microwave spectra at the decay phase 
are apparently better consistent with the kappa than 
Maxwellian fit, which does not allow to firmly con¬ 
clude that the Maxwellian distribution is maintained 
over the decay phase. One possible reason is that 
the actual energy distribution does deviate from the 
Maxwellian one at the decay phase due, perhaps, to 
nonthermal electron generation. It does not sound 
plausible, however, because this would imply that the 
particle acceleration did not take place at the stage 
of essential energy release, but appeared only at the 
decay phase, after the energy release is over. Instead, 
we believe it more likely that the plasma electron 
distribution remains (more or less) Maxwellian even 
at the decay phase; however, a nonuniform cooling 
results in a multi-temperature instead of the single¬ 
temperature plasma volume, so the deviation from 

^ The is noticeably smaller than one because estimates 
of the the statistical errors are overestimated as they are taken 
from the small antenna (less sensitive) but applied to big an¬ 
tenna (more sensitive) data. Here we are only interested in 
relative, rather than absolute, values of the measure. 


the Maxwellian fit at the decay phase is indicative of 
increasing source nonuniformity, rather than a non¬ 
thermal electron generation. This further implies 
that a uniform source with the kappa distribution can 
be used as a working model for nonuniform multi¬ 
temperature Maxwellian sources in other, more com¬ 
plex thermal solar flares. 

5. CONCLUSIONS 

We have described a nice solar flare whose spec¬ 
tral evolution in the microwave domain is consistent 
with a purely thermal behavior; only very minor de¬ 
viations (if any) from the Maxwellian distribution are 
consistent with the data. From this perspective this 
event is similar to a th ermal flare reported earlier by 
iGarv fc HnrfordI (jl989f ). whose spatial structures are, 
however, strikingly different from each other (a sin¬ 
gle source in the former case, while a complex source 
in our case). We believe that this difference comes 
from the different flare locations in these two cases. 
IGarv fc HnrfordI (jl989f ) reported a flare seen on disk, 
where the source of the gyro emission is observed at 
a quasi-parallel direction to the magnetic field, where 
no gyro-stripe structure is present. On the contrary, 
our event happened at the limb, so we likely observe 
a loop quasi-transversely to the magnetic field, when 
the spectral harmonic structure and the spatial gyro- 
stipes are the most pronounced. 

We detected a clear evolution of the source pa¬ 
rameters such as 5(t), T(t), and nth{t) suggesting 
a corresponding evolution of the magnetic and ther¬ 
mal energies at the radio source. However, the mea¬ 
sured energy components do not easily obey the en¬ 
ergy conservation requirement, implying that the ob¬ 
served source does not represent an isolated volume, 
but rather efficiently exchanges its energy with a big¬ 
ger surrounding volume. In particular, we concluded 
that the apparent decrease of the magnetic field at the 
radio source over the rise phase of the flare cannot be 
easily associated with the magnetic energy transfor¬ 
mation to the thermal or kinetic energy within a sin¬ 
gle isolated volume, but instead requires an upward 
propagating magnetic reconnection/plasma heating 
process, such as in the standard flare scenario with 
one remarkable difference, however: the absence of 
the nonthermal electron generation. We hope that 
the study of other thermal flares will better clarify the 
origin of such purely thermal events, which show a sig¬ 
nificant energy release observed through the plasma 
heating, but do not result in any measurable acceler¬ 
ation of the charged particles. 
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